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ABSTRACT 

We report on the first detection of extreme-ultraviolet (EUV) absorption variability in the 
Nevill AA770,780 mini-broad absorption line (mini-BAL) in the spectrum of the quasar 
(QSO) PG 1206H-459. The observed equivalent width (EW) of the Ne VIII doublet shows a 
^4cr variation over a timescale of 2.8 months in the QSO’s rest-frame. Both members of 
the Ne VIII doublet exhibit non-black saturation, indicating partial coverage of the continuum 
source. An increase in the Nevill covering fraction from fc = 0.59 ± 0.05 to 0.72 ± 0.03 
is observed over the same period. The Ne VIII profiles are too highly saturated to be suscep¬ 
tible to changes in the ionization state of the absorbing gas. In fact, we do not observe any 
significant variation in the EW and/or column density after correcting the spectra for partial 
coverage. We, thus, propose transverse motions of the absorbing gas as the cause of the ob¬ 
served variability. Using a simple model of a transiting cloud we estimate a transverse speed of 
1800 km s“^. Eor Keplerian motion, this corresponds to a distance between the absorber and 
the central engine of 1.3 pc, which places the absorber just outside the broad-line region. We 
further estimate a density of ^5x10® cm“^ and a kinetic luminosity of ~10^^-10"^^ erg s“^. 
Such large kinetic powers suggest that outflows detected via EUV lines are potentially major 
contributors to AGN feedback. 
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1 INTRODUCTION 

Broad absorption line (BAL; line spread. An ~fewx 1000 km s“^) 
and mini-BAL (An ~fewx 100 km s~^) systems detected via UV 
absorption troughs blue-shifted with respect to the QSO’s emis¬ 
sion redshift are thought to originate in an equatorial wind that 
is launched from the accretion disk near the central supermassive 
black hole (SMBH; e.g., Murray et al. 1995; Arav et al. 1995; Proga 
et al. 2000). The mass, momentum, and energy carried by such 
winds are increasingly invoked as the primary feedback mecha¬ 
nisms to explain the evolution of SMBH (Silk & Rees 1998; King 
2003, 2005) and their host-galaxies (Begelman & Nath 2005; Di 
Matteo et al. 2005). For example, it has been shown that an out¬ 
flow with a kinetic power of a few per cent of the Eddington lumi¬ 
nosity provides sufficient feedback to quench star formation, thus 
coupling the evolution of the host-galaxy and the growth of the 
central SMBH (Di Matteo et al. 2005). QSO outflows can also be 
responsible for metal enrichment in the interstellar medium (ISM), 
in the circumgalactic medium (COM), and even in the intergalac- 
tic medium (IGM; Wiersma et al. 2011). Therefore, the study of 
QSO outflows via BAL and mini-BAL systems is immensely im¬ 
portant. In particular, the incidence rate, outflow velocity, variabil¬ 
ity timescale, and covering fraction (fc) of the absorption troughs 


* E-mail: sowgatm@gmail.com 


provide essential clues about the location, geometry, and energetics 
of the flow. 

It is now well known that BAL systems show variability in 
absorption profiles over a time period of 1-10 yrs in the QSO rest- 
frame (e.g., Barlow 1994; Srianand & Petitjean 2001; Lundgren 
et al. 2007; Gibson et al. 2008; Capellupo et al. 2011; Filiz Ak 
et al. 2013; Vivek et al. 2012, 2014, and references therein). Filiz 
Ak et al. (2013) found a BAL fraction of ~50-60 per cent that 
exhibit variability over time-scales of 1-4 yrs in the QSO’s rest- 
frame. Additionally, BAL variability is found to be larger for longer 
timescales (e.g., Gibson et al. 2008; Capellupo et al. 2011; Vivek 
et al. 2014). For example, Capellupo et al. (2011) have found vari¬ 
ability in ~40 per cent of BAL absorbers on timescales of < 1 yr 
with that fraction increasing to ~65 per cent on timescales of 4-8 yr 
in the QSO rest-frame. BAL variability is generally explained by: 
(a) change in fc due to motions of the absorber across the line of 
sight and/or (b) change in the ionization state of the gas. In either 
case, it is possible to determine the dynamics and location of the ab¬ 
sorber with respect to the central engine under some assumptions. 
Note that the above explanations hold for mini-BAL variability as 
well. 

While BAL variability is well studied via several ongoing 
monitoring campaigns by different groups with a fairly large num¬ 
ber of objects, very little is known for mini-BALs. In fact, variabil¬ 
ity in mini-BAL absorption is studied only in a handful of cases 
(i.e., Hamann et al. 1997; Narayanan et al. 2004; Misawa et al. 
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2007; Rodriguez Hidalgo et al. 2011; Misawa et al. 2014). Inter¬ 
estingly, Misawa et al. (2007) have demonstrated that the afore¬ 
mentioned causes for line variability cannot reproduce the rapid, 
simultaneous variation in multiple troughs of the mini-BAL sys¬ 
tem detected towards HS 1603-1-3820. Instead, they proposed that 
a clumpy “screen” of gas with variable optical depth and covering 
fraction, located between the continuum source and the absorbing 
gas, mimics the variation in the ionizing continuum. The variability 
study of a sample of seven mini-BALs was recently presented by 
Misawa et al. (2014). Assuming that the observed variations result 
from changes in the ionization state of the mini-BAL gas, they have 
inferred lower limits to the gas density of ~10® — 10® cm~®, corre¬ 
sponding to upper limits on the distance of the absorbers from the 
central engine of ~ few kpc. Note that unlike BAL systems, absorp¬ 
tion line parameters can be more robustly determined for mini-BAL 
systems because (a) the QSO continuum is easily identifiable even 
in the presence of mini-BALs, (b) doublets are not self-blended due 
to relatively narrow line profiles, and (c) the kinematic structure 
is often discernible due to lack of heavy line-saturation (see e.g., 
Muzahid et al. 2012). Mini-BAL systems are thus excellent tools 
for measuring outflow properties, e.g., density, metallicity, outflow 
mass, mass-flow rate, and kinetic luminosity that are essential for 
quantifying the so called “AGN feedback” (Arav et al. 2013). 

Recently, in an archival Hubble Space Telescope (HST)/ 
Cosmic Origins Spectrograph (COS) study, we have reported a 
new population of associated absorbers detected via absorption 
from EUV ions (e.g., NeVIII, NalX, and MgX) in the spectra 
of intermediate-^ QSOs (Muzahid et al. 2013, hereafter M13). 
We have found that a significant fraction (~40 per cent) of the 
intermediate-^: QSOs exhibit outflows observable via EUV lines. 
These absorbers are very highly ionized with ionization parameters 
(i.e. 0.5< log U <1.0) consistent with their being intermediate to 
X-ray “warm absorbers” (Halpern 1984) and BALs. The velocity 
spreads (100-800 km s“^) of the NeVIII absorption lines suggest 
that they are mini-BALs. Here we present a fortuitous case from the 
sample presented in M13 in which the QSO (PG 1206-1-459) was 
observed at three different epochs. PG 1206-1-459 is an UV-bright 
QSO with a U-band magnitude of ~15.4 and has an emission red- 
shift of Zsni ~ 1.163. The multi-epoch spectroscopic observations 
of the QSO allowed us to study variability of the OIV, NIV, and 
Ne VIII absorption. We found a ~4cr variation in the Ne VIII equiv¬ 
alent width. This is the first case of EUV mini-BAL variability re¬ 
ported so far. This article is organized as follows: In Section 2 we 
present observations and data reduction of the QSO PG 1206-1-459. 
Absorption line analysis and the QSO’s properties are presented in 
Section 3. In Section 4 we discuss the implications of our findings. 
In Section 5 we briefly summarized our main results. Throughout 
this article we adopt a flat ACDM cosmology with Ha = 0.73, 
Om = 0.27, and Hq = 71 km s“^ Mpc“^. Unless specified 
otherwise, the EWs are given in the observed frame whereas the 
timescales are in the QSO’s rest-frame. 


2 OBSERVATIONS AND DATA REDUCTION 

Far-UV (EUV) spectra of the UV-bright QSO PG 1206-H459 (Zem 
= 1.163) were obtained using HST/COS Cycle-17 observations, 
under program ID: GO-11741 (PI: Todd Tripp). These observations 
consist of G130M and G160M grating integrations at a medium 
spectral resolution of R ~20,000 (FWHM ~18 km s“^). The data 
were retrieved from the HST archive and reduced using the STScI 
CALCOS v2.21 pipeline software. The reduced data were flux cal- 


Table 1. Details of the HST/COS observations. 


Grating"' 

(A) 

Tilt“ 

‘'exp 

(sec) 

Date of 

Observation 

MJD“ 

Epoch S/N^ 

G160M 

1600 

1 

2361.2 

2009-12-29 

55194.5 

1 


G160M 

1600 

2 

3137.2 

2009-12-29 

55194.5 

1 


G160M 

1623 

3 

3137.2 

2009-12-29 

55194.5 

1 


G160M 

1623 

4 

3137.2 

2009-12-29 

55194.5 

1 

11 

G160M 

1600 

1 

3137.2 

2010-01-05 

55201.5 

2 


G160M 

1600 

2 

3137.2 

2010-01-05 

55201.5 

2 


G160M 

1600 

3 

3137.2 

2010-01-05 

55201.5 

2 


G160M1 

1600 

4 

3137.0 

2010-01-05 

55201.5 

2 

10 

G160M 

1623 

1 

2361.2 

2010-06-14 

55361.5 

3 


G160M 

1623 

2 

3137.1 

2010-06-14 

55361.5 

3 


G160M 

1623 

3 

3137.2 

2010-06-14 

55361.5 

3 


G160M 

1623 

4 

3137.2 

2010-06-14 

55361.5 

3 

12 

G130M 

1309 

1 

2406.2 

2010-01-04 

55200.5 

2 


G130M 

1309 

3 

3137.2 

2010-01-04 

55200.5 

2 


G130M 

1309 

4 

3137.2 

2010-01-04 

55200.5 

2 


G130M 

1327 

1 

3137.2 

2010-01-04 

55200.5 

2 


G130M 

1327 

3 

3137.2 

2010-01-04 

55200.5 

2 


G130M 

1327 

4 

2406.2 

2010-01-05 

55201.5 

2 

27 


Notes- “Grating used. ^’Central wavelength. “Fixed-pattern offset position. 
'^Exposure time in second. “Modified Julian date corresponding to the 
observation date. 7 S'/A per resolution element of the coadded data at 
~1575 A for G160M grating and at ~1230 A for G130M grating. ^This 
exposure {IbloOSkqq) was bad and not used for co-addition. 

ibrated. To increase the spectral signal-to-noise ratio (S/N), spec¬ 
tra from individual exposures were aligned and coadded using the 
IDL code “coadd_xld” developed by Danforth et al. (2010). The 
properties of COS and its in-flight operations are discussed by Os- 
terman et al. (2011) and Green et al. (2012). The coadded spectra 
are significantly oversampled, with six raw pixels per resolution el¬ 
ement. We thus binned the data by six pixels. This improves the 
spectral S/N per pixel by a factor of y/6. All our measurements 
and analyses were performed on the binned data. The binning does 
not affect our EW and/or fc measurements since the absorption 
lines of interest are much broader than the spectral resolution. Con¬ 
tinuum normalization was done by fitting the line-free regions with 
smooth low-order polynomials. Since there is no Lya forest crowd¬ 
ing at low-a and the absorption lines of interest do not fall on top 
of any emission lines, the continua were very well defined in most 
part. COS wavelength calibration is known to have uncertainties 
at the level of 10-15 km s“^ (Savage et al. 2011; Meiring et al. 
2013). Using molecular hydrogen absorption lines, Muzahid et al. 
(2015b) have noted that while in majority of the cases wavelengths 
are accurate within ~ ±5 km s“^, large offsets (>10 km s“^) are 
often noticeable at the edges of the COS spectra. We noticed that 
the wavelength calibration uncertainty is usually <5 km s“^ for 
PG 1206-1-459. However, the uncertainty increases to ~10 km s~^ 
at the edge of the spectra as also noted by Muzahid et al. (2015b). 

Besides COS, we have used the HST/Space Telescope Imag¬ 
ing Spectrograph (STIS) spectrum of PG 1206-1-459. The STIS 
spectrum was obtained using the E230M grating under program 
ID: 8672 (PI: Chris Churchill). We refer the reader to Ding et al. 
(2003) for the full details on the STIS observations. The spectrum 
is only used for calculating rest-frame 1350 A flux. 

In Table 1 we summarize the details of the observations. It 
is clear from the table that the G160M observations were taken 
in three different epochs. Epoch-1 and Epoch-2 are separated by 
~7 days, corresponding to ~3.2 days in the QSO’s rest-frame. 
Epoch-1 and Epoch-3 are separated by ~6 months (~2.8 months 
in the QSO’s rest-frame). We note that the G160M spectrum cov- 
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-19.6 -19.2 -18.8 


Outflow Velocity (10^ km s“^) 


Figure 1. Normalised flux against the outflow velocity, measured with respect to the Zem = 1.163, for different ions. The absorption from the outflow are 
shown by the shaded regions. All other absorption are unrelated and intervening. The blended pixels are marked by the horizontal red bars. In each panel the 
data with different colors represent different epochs as indicated in the top. No significant variation is seen in the left panel between Epoch-1 and Epoch-2 for 
any of the troughs. We therefore co-add the Epoch-1 and Epoch-2 data to obtain a better S/N spectrum which is represented by Epoch-1' in the right panel. 
While no appreciable valuation is seen between Epoch-1' and Epoch-3 in the N IV and O IV absorption, the Ne VIII doublet shows a significant change in the 
absorption troughs. Note that only one central wavelength setting was used for the GI60M grating obseiwations in Epoch-2 and the 0 IV line falls in the gap 
of the Epoch-2 spectrum. 


ers only the Niv A765, OlV A 787, and NeVIII AA770,780 tran¬ 
sitions from the mini-BAL outflow we studied here. Unfortunately, 
other important high-ionization species (e.g., NalX AA681,694, 
MgX AA609,624) fall in the wavelength range covered hy the 
G130M grating and are observed only once. 

3 ANALYSIS 

3.1 Properties of the absorber 

The mini-BAL outflow from QSO PG 1206-1-459 (Zem = 1.163) 
was first reported in Ml3 as a part of our HST/COS sur¬ 


vey of intrinsic absorbers detected via Ne VIII AA770,780 dou¬ 
blet transitions. The absorption redshift of z^bs = 1.02854 cor¬ 
responds to an outflow velocity of ~19,200 km s“^. Both mem¬ 
bers of the NeVIII doublet are detected with EW/(tbw > 20. 
The Ne VIII absorption is spread over Av ~360 km s~^ with a to¬ 
tal column density of N{Ne VIIl) > 10^® cm“^. A very weak Lya 
line is detected with log(Af(Hl)/cm“^) <14 in the HST/STIS 
spectrum. Interestingly, no singly or doubly ionized species are de¬ 
tected from this absorber. Besides NeVIII, we have reported the 
presence of several other high-(Ar VIII, NalX, MgX) and low- 
(OIV, NIV, O V, N V) ionization species (see Fig. 7 of M13). But 





































































































4 S. Muzahid et al. 


only N IV A765, Ne VIII A770,780, and O IV A787 lines have multi¬ 
epoch coverage. 

In Fig. 1 we show the absorption profiles of NIV A765, 
OlV A 787, and NeVIII AA770,780 transitions for the three dif¬ 
ferent epochs. No significant variation in any of these absorption 
troughs between Epoch-1 and Epoch-2 (separated by ~3.2 days in 
the QSO’s rest-frame) is apparent from the figure. To quantify the 
degree of variability in an absorption, we measure EWs at different 
epochs and define absorption variability significance as follows: 


Sij = 


EWj - EWi 




+ ) 


( 1 ) 


where, EWj and crBWj are the equivalent width and its error, re¬ 
spectively, as measured in j-th epoch. Note that a positive 5'ij indi¬ 
cates an increase in EW at later epoch. We have estimated that the 
error in the observed EW due to the continuum placement uncer¬ 
tainty is on the order of ctew. We, thus, introduce the v/2 factor in 
the denominator in Eq. 1 to take care of the continuum placement 
uncertainty (see also Misawa et al. 2014). 

Next, we choose a set of 39 strong absorption lines that are 
detected with EW/(Jew > 5 in all three epochs in the wavelength 
ranges 1440-1587 A and 1623-1775 A^. Details of these lines are 
given in Appendix A. Note that all but three of these lines, i.e., 
NIV A765, NeVIII A770 and A780 lines from the mini-BAL^, are 
intervening and should not show any considerable variation with 
time. For all these absorption lines we have calculated S'ij. The top 
panel of Fig. 2 shows the distribution of 5'ij as measured between 
Epoch-1 and Epoch-2 ( 5 i 2 ). It is apparent from the figure that none 
of the absorption lines, including the three from the mini-BAL sys¬ 
tem, shows variability with a significance of >2a. A Gaussian fit to 
the 5 i 2 distribution gives a mean of fj, = —0.3 ±0.2 and a standard 
deviation of ct = 1.0 ± 0.2. A 2a variation, thus, is not significant 
for the data we have. Since there is no considerable variation in any 
absorption we co-add the Epoch-1 and Epoch-2 data in order to 
increase the spectral S/N. The resultant data, with marginally im¬ 
proved S/N of ~14 per resolution element, are shown as Epoch-1' 
in Eig. 1. The Epoch-1' spectrum has spectral coverage of 1150- 
1800 A. Eor interested readers we have shown the Epoch-1' spec¬ 
trum in Appendix B. 

In the bottom panel of Fig. 2 we show the 5ij distribution as 
calculated between Epoch-1' and Epoch-3 (Si' 3 ). Here also we 
find that none, but two of the absorption lines, shows variability 
with a significance of >2a. A Gaussian fit to the S 113 distribution 
results in similar fj, and a values as we obtained for the 5 i 2 dis¬ 
tribution. The two absorption lines showing ~4a variation are the 
Ne VIII A770 and A780 of the mini-BAL outflow (see Appendix A). 
Clearly, the variations are statistically significant and cannot be at¬ 
tributed to some systematic uncertainties. The N IV line, however, 
does not show any considerable variability between Epoch-1' and 
Epoch-3, e.g., shows only a —OAa variation. Moreover, the OlV 
line shows only a ~0.2a variation between Epoch-1 and Epoch-3. 

In a situation when an absorber covers the background contin¬ 
uum emitting region only partially with a covering fraction of fc, 
then the observed and “true” optical depths (ta) are related via: 

exp (-Ta'"'’) = = (1 - /c) ± fc exp{-Tx'^^) . ( 2 ) 

Here is the observed normalized flux. Here we assume that 


^ We do not use lines within ~20 A from both blue- and red-end of the 
spectra due to wavelength calibration uncertainty. 

^ O IV A787 is not considered here since it is not covered in Epoch-2 data. 



Figure 2. The 5'ij distributions between Epoch-1 & Epoch-2 {S 12 , top) and 
between Epoch-1' & Epoch-3 (5i/3, bottom) for a set of 39 strong absorp¬ 
tion lines detected in all three epochs (see text). The details of these absorp¬ 
tion lines are summarized in Appendix A. Each of the observed 5ij distri¬ 
butions is modelled as Gaussian with fit parameters given in each panel. No 
absorption lines shows a >2a variability between any of the epochs except 
for the NeVIII AA770,780 transitions which show ~4c7 variability, in the 
bottom panel, between Epoch-1' & Epoch-3. 


the background source is partially covered by a homogeneous ab¬ 
sorber with a constant fc. However, fc can be strongly dependent 
on velocity along the absorption trough (e.g., Arav et al. 2008). In 
the case of a doublet with rest-frame wavelengths of Ai and A 2 , and 
oscillator strengths of fi and / 2 , the normalized fluxes of and 
at any velocity with respect to the line centroid are related by: 


Rt^iv) = (1 - fc) ± /c X 


/- Rf^{v)-l + fc 

V fc 


7 


(3) 


where, 7 = / 2 A 2 // 1 A 1 , which is close to 2 for the NeVIII dou¬ 
blet (see e.g., Petitjean & Srianand 1999; Srianand & Shankara- 
narayanan 1999). Using Eq. 3 we have calculated the NeVIII cov¬ 
ering fraction as a function of outflow velocity as shown in Fig. 3 
(see also Muzahid et al. 2012, 2013). The median fc value obtained 
for the Epoch-1' data is 0.59 ± 0.05, whereas it is 0.72 ± 0.03 for 
the Epoch-3 data. The error bars correspond to the standard devi¬ 
ations of fc values measured in each pixel. Note that the median 
fc values are determined from the core of the absorption, i.e., only 
using the pixels with highest optical depths marked by the horizon¬ 
tal bars. Nevertheless, it is evident from the figure that the Epoch-3 
data show higher fc values as compared to the Epoch-1' data in 
almost every pixel over the entire velocity range. 

Einally, in order to evaluate whether the Ne VIII column den¬ 
sity has increased with time, we perform EW measurements in the 
covering fraction corrected spectra (i.e. the spectra obtained by in¬ 
verting the Eq. 2). For both the Epoch-1' and Epoch-3 data we mea¬ 
sure EWobs(NeVIII A770) = 2.00 ± 0.04 A. Thus, no significant 
change in Ne VIII equivalent width is noticeable in the partial cover¬ 
age corrected spectra. A change in N(Ne VIIl) without a significant 
change in EW is possible in the flat part of the curve-of-growth. 
















Covering Fraction 
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Outflow Velocity (km s“*) 

Figure 3. The Ne VIIl covering fraction against the outflow velocity (rela¬ 
tive to the Zem) for the two different epochs. The thick blue curve represents 
Epoch-1', whereas the thin red curve is for Epoch-3. The median /c values 
as shown by the horizontal dashed lines exhibit an increase from 0.59±0.05 
(std) to 0.72 ± 0.03 (std) from the Epoch-1' to Epoch-3. The height of the 
shaded horizontal bars represent the Icr range in the measured /c values. 
Note that the fc measurements in the grey shaded region are not reliable 
due to an unknown blend in the red wing of the Ne VIlI A780 absorption. 

However the fact that we do not see any change in the wings of 
the profile indicates no considerable change in A'^(Ne VIIl) over Ihe 
course of 2.8 months in the QSO’s rest-frame. 

3.2 Properties of the background QSO 

The UV bright QSO PG 1206-1-459 has a l/-band magnitude of 
~15.4. The Catalina Real-Time Transient Survey (CRTS; Drake 
et al. 2009) light-curve, shown in Fig. 4, indicates very little fluc¬ 
tuation in the H-band magnitude and hence in the continuum flux. 
In addition, no significant variation in the UV continuum is seen 
in the COS data. We found a maximum of only ~30 per cent vari¬ 
ation in the UV flux between Epoch-1' and Epoch-3. Using the 
Mg II emission line width, Chand et al. (2010) have estimated a 
black hole mass of Mbh ~1O®M0 for this QSO which corre¬ 
sponds to a Schwarzschild radius of 7?sch ~3 x 10^^ cm. Fol¬ 
lowing Hall et al. (2011), we estimate the diameter of the disk 
within which 90 per cent of the 2700 A continuum is emitted, 
D 2700 ~50i?sch ~1.3 X 10^® cm®. From the rest-frame 1350 A 
flux (i.e., ~5xl0“®® erg cm“^ s~® A“®), measured from the 
STIS spectrum, we obtain the broad line region (BLR) size of 
Rbl,r ~0.7-0.9 pc by adopting the R — L(1350A) relation given 
in Kilerci Eser et al. (2015). 


4 DISCUSSION 

We report on the variability in Ne VIII absorption originating from 
a mini-BAL outflow from the QSO PG 1206-1-459 (Zem = 1.163). 
The absorber is detected at Zabs = 1.02854 corresponding to an 
outflow velocity of ~19,200 km s“®. The mini-BAL system, with 
NeVIII absorption kinematically spread over Av ~360 km s~® 
and with a total column density of log (W(Ne VIIl)/cm~^) > 16, 
was first reported in M13 (see their Table 2). Here we report a ~4 (t 

® We note that the Mbh and the rest-frame 3000 A luminosity of the 
present source are, respectively, ~3 times lower and ~5 times higher than 
the one studied by Hall et al. (201 1). The facts ensure that the disk temper¬ 
atures, and hence the size of the continuum emitting regions, are consistent 
with each other at ~15 per cent level. 


variation in the observed EWs of the NeVIII A770 and A780 ab¬ 
sorption lines over a timescale of 2.8 months in the QSO’s rest- 
frame. Such variations cannot be attributed to systematic uncer¬ 
tainties since the neighboring unrelated, intervening absorption do 
not vary with such a high significance. This is the first reported 
case of an EUV mini-BAL variability. However, we did not ob¬ 
serve any considerable variation in the Ne VIII absorption troughs 
on a timescale of ~3.2 days in the QSO’s rest-frame. Additionally, 
two other low-ionization species (OIV and NIV) with multi-epoch 
observations do not exhibit any appreciable variation at any time, 
e.g., show < 2a variations. Note that the mini-BAL absorbers in 
the sample of Misawa et al. (2014) show < 2cr variability over a 
similar timescale. Thus a ~4(j absorption variability in Ne VIII over 
2.8 months in the QSO’s rest-frame is somewhat extreme. How¬ 
ever, we point out that the mini-BAL absorbers in their sample are 
detected via Si IV, CIV, and N V absorption lines which have ion¬ 
ization potentials similar to NIV and OIV. The lack of significant 
variability in NIV and OIV lines in the present system, therefore, 
is consistent with the observations of Misawa et al. (2014). 

As mentioned in Section 1, the change in ionization state of 
the absorber and the change in fc due to absorber’s motions are 
the two most plausible reasons for absorption variability. Below we 
discuss the two scenarios in detail. 


4.1 Variability due to change in ionization 

Variation in the optical depth of an absorption line (i.e. variability) 
naturally occurs due to the change in ionization state of the absorb¬ 
ing gas. The change in ionization state of photoionized gas can be 
induced via a change in the ionizing continuum. In the “screening 
gas” scenario of Misawa et al. (2007), the change in optical depth 
of “screening gas” can also lead to absorption line variability. Nev¬ 
ertheless, variation in the continuum flux is thought to be the most 
probable reason for the variation in the optical depth of “screening 
gas”. Here we have found that the U-band magnitude and/or the 
UV continuum of the background QSO do not show any signifi¬ 
cant variation (Section 3.2). The weak variation (up to 30%) seen 
in the UV continuum will not change the Ne VIII column density 
appreciably. This is because A^(Ne VIIl) reaches at its peak around 
ionization parameter of log U ~ 1 and shows a very little change 
over a wide range in log U around the peak, i.e., log U ~0.5-L5 
(see Fig. 12 of M13 for example). Thus a maximum of 30% change 
in log U induced by the change in the UV continuum, assuming 
the shape of the ionizing spectrum is not changing, will not be dis¬ 
cernible. Furthermore, the total column density of Ne VIII as given 
in M13 (i.e. log {N/cm~^) > 16, after correcting for the partial 
coverage) indicates that the profiles are too heavily saturated to be 
susceptible to the change in the ionization state. Finally, the lack 
of any considerable variation in the “true” NeVIII EWs, as mea¬ 
sured from the partial coverage corrected spectra, and the lack of 
any change in the wings of the NeVIII profile (Section 3.1) indi¬ 
cates that it is highly unlikely that the observed line variation is due 
to the change in ionization state of the absorbing gas. 

4.2 Variability due to absorber’s motion 

Motions of the absorbing gas cloud, transverse to the line of sight, 
could result in a change in fc and hence in the “apparent” optical 
depth of absorption. In Section 3.1 we have shown that the fc of the 
gas phase giving rise to the Ne VIII absorption shows an increase 
from fc = 0.59±0.05 to 0.72±0.03 over a timescale of 2.8 months 
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Figure 4. The CRTS V-band magnitudes of QSO PG 1206+459 plotted 
against the modified Julian date (MJD). Clearly, the background source 
shows a negligible (< 0.1 mag) variation over the entire MJD range. The 
Epoch-1 and Epoch-3 are indicated by the arrows. 

in the QSO’s rest-frame. Such an increase in /c value can occur 
primarily because of the two reasons: (i) transverse expansion of 
the absorbing gas cloud, and/or (ii) absorber’s bulk motions across 
the line of sight. 

4.2.1 Expansion across the line of sight 

If the entire absorber is already within our cylinder of sight then 
any transverse motion will lead to a decrease in /c in contrast to 
what we observe here. Nevertheless, a gas cloud, without a sig¬ 
nificant transverse motion, can expand in size and consequently 
increase the covering fraction. In this scenario, in Epoch-1, the 
transverse size of the absorbing gas cloud is estimated to be rj_ = 
v/7c-^270o/2 ~5.0 X 10^® cm. An increase in radius from ri to r 2 , 
due to expansion, will increase the covering fraction from /i to /2 
via the relation: ri/r 2 = sjf\/f 2 . The fact that the Ne VIII cov¬ 
ering fraction increases from ~0.59 to ~0.72 in the rest-frame 2.8 
months allows us to estimate an expansion rate of ~720 km s“^. If 
this were adiabatic expansion at the speed of sound then it would 
imply a temperature of ~4 x 10^ K. But such a high temperature 
is incompatible with the ionization states of metals that we observe 
in the absorber. 

More than 14 absorption troughs from eight different low- 
(e.g. O IV, N IV, N V) and high-(e.g. Ne VIII, Na IX, Mg X) ions are 
detected from this system (see Fig. 7 of Ml 3). Interestingly, in Ml 3 
we have shown that these ions exhibit ionization potential (IP) de¬ 
pendent covering fractions with higher ions covering more of the 
continuum emitting region. For example, the N IV and O IV have 
/c values of ~0.4 and 0.2, respectively, as opposed to ~0.6 we 
measured for the Ne VIII. Such observations clearly suggest mul¬ 
tiphase structure of the absorber with different ions having differ¬ 
ent projected area across the line of sight. Different phases of a 
multiphase photoionized absorber are often found to have very dif¬ 
ferent densities and temperatures causing pressure gradients in the 
absorbing gas (e.g. Muzahid et al. 2012). Such pressure gradients 
could, eventually, lead to an expansion. 

It is evident from Fig. 12 of M13 that A'(NeVIIl) peaks 
at log (7 = 1.0, whereas, N{Ol'V) and A’(Niv) peaks at 

log U = —1.0. Note that the photoionization equilibrium temper¬ 
atures for the high- and low-ionization phases are T ~10® and 
~10^ K, respectively. At these temperatures, the gas-pressures for 
both the phases are at least 2 orders of magnitude higher than the 
corresponding radiation pressures, for the whole range of densi¬ 
ties. Therefore, the total pressure, dominated by the gas-pressure. 


of the low-ionization phase is a factor of ~10 higher than that of 
the high-ionization phase. As a consequence, it is natural for the 
low-ionization gas phase to exhibit a more rapid expansion than 
the NeVIII bearing gas. However, we do not find any appreciable 
variability in the NIV and/or OlV absorption. Hence, the expan¬ 
sion scenario, although highly plausible in case of a multiphase 
absorber, cannot explain our observations. In view of the above ar¬ 
guments, we disfavor the expansion interpretation for the change in 
/c in this system. 

4.2.2 Bulk motions across the line of sight 

In a scenario in which an absorber is entering to our cylinder of 
sight, due to a rotational motion around the central SMBH, can ex¬ 
plain the phenomenon of increase in /c with time. In Fig. 5 we 
present a simple model in which a cloud with a radius of Rc is 
transiting a continuum emitting disk with a radius of Rd. Since the 
absorber covers ~60-70% of the background source, we reason¬ 
ably assume Rc ~ Rd- The common area between the cloud and 
the disk, under such an approximation, can be expressed as 

A = 2Rjcos-\d/2Rd)-^ds/(2Rdf^, (4) 
leading to a covering fraction of 

fc = - cos“^(«;/2) - - s/1 - «:2/4 . (5) 

TT TT 

Here, d is the center-to-center projected separation between the 
disk and the absorber, and k is a dimensionless parameter defined 
as K = d/Rd. Note that d can vary from 0 to 2Rd and k can 
have values between 0 and 2. The right panel of Fig. 5 shows 
how /c changes with n. The observed increase in fc value from 
0.59 to 0.72, as shown by the horizontal dotted lines, corresponds 
to Ak ~0.2 leading to Ad ~Q.2Rd = O. 2 D 2700/2 ~1.3 x 
10^® cm. The transverse speed of the absorber is therefore, vj_ — 
Ad/At ~1800 km s“^. Note that the transverse speed is ~10 
times smaller than the radial outflow velocity we measured from 
the spectra (i.e. ~19,200 km s”*^). If we assume that the vj_ can 
be matched to a Keplerian velocity (see e.g., McGraw et al. 2015), 
then the derived distance between the absorber and the central en¬ 
gine is rd ~4.0x 10^® cm ~ 1.3 pc. This suggest that the outflowing 
gas is located just outside the BLR region. With the observed radial 
velocity the cloud will move only Ar^ ~4.5 x 10“® pc in the ra¬ 
dial direction over the course of the 2.8 months. Clearly, the Avd is 
too small, as compared to the rd, to cause any significant change in 
the photon flux due to the change in the distance between the QSO 
and the absorber. 

A detailed photoionization model for this absorber was pre¬ 
sented in M13. Here we recall the important model parameters: 
i.e., log U ~1, metallicity of [X/H] ~1.0, and total hydrogen col¬ 
umn density of log (Wn/cm”^) = 20.7. Inserting rd = 1.3 pc 
in Eq. 4 of M13^, we calculate the density of the NeVIII bear¬ 
ing gas to be nn ~5.3 x 10® cm“®. Given the density, the line 
of sight thickness of the cloud turns out to be Sr = Wh/tih 
~ 9.5 X 10^® cm. On the other hand, since the absorber covers 
at least ~60 per cent of the continuum emitting region, the mini¬ 
mum transverse size of the absorber is r± — s/J^D 27 oo 12 ~5.0 x 

( lowO = ^“2 A ■ (t^) ' - (T - 1.25, 

where log o = 32.1, is the monochromatic luminosity of the OSO at 
912 A 

the Lyman continuum in units of 10^® erg s“^ Hz“^. 
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Figure 5. Left: A simple model of a transiting cloud. The cloud with a radius Rc is crossing the disk with a radius R^. The cloud’s center is separated by a 
distance of d from the center of the disk. The covering fraction of the cloud is the ratio of the area common to both the circles (A) to the total area of the disk. 
Right: The covering fraction as a function of re = djR^, assuming Rc ~ Rd- 


10^® cm. Evidently, r± > 50 Sr, which is suggestive of a “sheet¬ 
like” geometry. 

In a standard thin-shell model, the kinetic luminosity of an 
outflow can be expressed as Ek = 47r/impCnC/A^HfdU^, where, 
trip is the proton mass, r. = 1.4 accounts for the mass of he¬ 
lium, and Cn and Cf are, respectively, the global and local cov¬ 
ering fractions (see Appendix C of Muzahid et al. 2015a). In M13 
we have shown that 40 per cent of the intermediate redshift QSOs 
exhibit highly ionized outflows detected via the NeVIII absorp¬ 
tion. Therefore, Cn can be taken as 0.4 for the general popula¬ 
tion. However, we note that the outflow velocity (~19,200 km s~^) 
we measure for this system is very rare. Only one out of 20 
QSOs shows such a high outflow velocity. We, thus, conservatively 
use Cn = 1/20 for our calculations. The Cf is related to the 
dumpiness of the absorbing material. For simplicity, Cf is taken 
to be unity for the highly ionized diffuse gas. Using the distance 
Td ~1.3 pc, derived from our simple model, we estimate a kinetic 
luminosity of Ek ~2.1 x 10"^® erg s“^. The Ek value increases 
to ~1.7 X 10"'“* erg s“^ if we assume Cn = 0.4. Note that the 
Ek of erg s~^ is among the highest kinetic luminosities 

measured till date (see Table 10 of Dunn et al. 2010). Such large ki¬ 
netic luminosities suggest that the outflows detected via the Ne VIII 
absorption are potentially major contributors to AGN feedback. 

It is now important to address whether we could explain the 
lack of variability in the low-ionization lines, i.e. OIV and N IV, 
with our simple model. Naively one would expect that the low- 
ionization gas phase would trace a high density, compact region. 
In M13, using OlV A608 and A787 transitions, we determined the 
OIV covering fraction to be ~0.2. Therefore, the transverse size of 
the low-ionization gas phase, giving rise to OIV (and N IV) absorp¬ 
tion, is only a factor of ~2 smaller than the Ne VIII bearing gas. 
The fact that all the high- and low-ionization lines coincide in ve¬ 
locity, possibly, indicates that the gas phases are co-spatial with the 
low-ionization lines stemming from the core of the absorber with a 
lower fc value. While such a scenario is reasonable, the difference 
in transverse sizes is not quite significant. 

We recall that according to Eq. 1, the change in EW, and hence 
in Sij, depends on the EW itself. From Fig. 1 it is evident that the 


both OIV (EW= 0.37 ±0.05 A) and Niv (EW = 0.74± 0.03 A) 
lines are much weaker than the NeVIII (EW = 1.18 ± 0.04 A). 
Thus, one possibility here is that variations in those weak lines, 
OIV in particular, are below our detection limit. In order to check 
that, we estimated the S^j value for OlV in the partial coverage 
corrected spectra. However, we did not find any significant varia¬ 
tion. Higher S /N multi-epoch observations are essential for reach¬ 
ing a firm conclusion on the degree of variability of these ions. As 
noted earlier, the lack of appreciable variation in the low-ionization 
species over a timescale of 2.8 months, however, is consistent with 
the recent observations by Misawa et al. (2014). 


5 SUMMARY 

We have studied the variability of the NeVIII AA770, 780 dou¬ 
blet originating from a mini-BAL outflow with an ejection veloc¬ 
ity of ~19,200 km s"i from the UV bright QSO PG 1206-1-459 
(2em = 1.163). This is the first case of mini-BAL variability de¬ 
tected in EUV transitions. The NeVIII doublet shows variability 
with a ~4 ct significance over a timescale of 2.8 months in the 
QSO’s rest-frame. However, no significant variation is observed 
over a shorter timescale of 3.2 days. Additionally, we did not ob¬ 
serve any appreciable variation in lower ionization O IV and NIV 
lines. 

We have explored three different possibilities for the observed 
variability: (a) change in the ionization conditions, (b) expansion 
of the absorber, and (c) bulk motion of the absorbing gas across 
the line of sight. From the lack of variability in the partial cover¬ 
age corrected spectra and estimated large column density of Ne VIII 
(N > 10^® cm~^) we rule out the possibility of a change in ioniza¬ 
tion conditions. The expansion scenario leads to an unreasonably 
high gas temperature that is inconsistent with the ionization state 
of the detected metal lines. Besides, we note that the low-ionization 
lines should exhibit a more rapid expansion than the Ne VIII in con¬ 
trast to what we observed here. We favored a scenario in which the 
absorbing gas is entering into the cylinder of sight and hence show¬ 
ing an increase in covering fraction. Using a simple model of a 
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transiting cloud, with a size on the order of the continuum emitting 
region, we derived the transverse speed, distance from the central 
engine, density, and kinetic luminosity of the outflow. In order to 
explain the lack of variability in the low-ionization lines, we pre¬ 
ferred a multiphase structure of the absorbing gas in which the low- 
ionization species stem from the core of the absorber with lower /c 
values and presumably already in the cylinder of sight. Nonethe¬ 
less, the transverse size of the low-ionization gas phase is only a 
factor of ~2 smaller than that of the high-ionization phase, posing 
an apparent tension to our favored scenario. High S/N, multi-epoch 
observations are needed for a firm conclusion on the variability of 
the weak, low-ionization lines. 

Eight different ions with a wide range in IPs (i.e. 40-400 eV) 
are detected in this absorber (see Fig. 7 & 8 of M13). A study of 
both short-term and long-term variability with multi-epoch HST 
observations, in future, for more than 10 absorption troughs would 
provide further insights on the dynamics and the structure of this 
intriguing mini-BAL outflow. 
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Figure Bl. Epoch-1' spectrum of PG 1206+459 showing full wavelength coverage obtained using COS FUV (G130M+G160M) gratings. Error in each pixel 
is shown in blue. Continuum is shown by the smooth red curve. The absorption lines studied here are also marked. The wavelength ranges affected by the 
geo-coronal emission lines are indicated by the horizontal bars. 
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Table Al. Summary of EWs and S'ij measurements. 


Ai ( A ) 

(1) 

^2 ( A ) 

(2) 

EWi ( A ) 
(3) 

EW 2 ( A ) 

(4) 

5' i 2 

(5) 

EWi / ( A ) 
(6) 

EWg ( A ) 

(7) 

5 i '3 

(8) 

1460.7441 

1462.1398 

0.23±0.02 

0.21±0.02 

-0.6 

0.20±0.02 

0.16±0.02 

-1.3 

1463.0298 

1464.7288 

0.44±0.02 

0.43±0.02 

-0.5 

0.45±0.02 

0.42±0.02 

-0.9 

1464.6884 

1467.8639 

1.27±0.03 

1.26±0.03 

-0.1 

1.24±0.02 

1.24±0.03 

0.0 

1468.5472 

1470.6407 

0.85±0.03 

0.79±0.03 

-1.7 

0.82±0.02 

0.77±0.02 

-1.7 

1471.0566 

1472.1326 

0.62±0.01 

0.59±0.02 

-1.4 

0.64±0.01 

0.62±0.01 

-1.2 

1474.4214 

1475.6957 

0.91±0.01 

0.89±0.02 

-1.1 

0.91±0.01 

0.91±0.01 

-0.5 

1475.6957 

1477.9813 

1.17±0.02 

1.17±0.03 

0.0 

1.18±0.02 

1.19±0.02 

0.4 

1479.3688 

1480.5581 

0.13±0.02 

0.19±0.02 

1.8 

0.13±0.02 

0.15±0.02 

0.7 

1481.4803 

1483.7620 

0.43±0.03 

0.44±0.03 

O.I 

0.45±0.02 

0.44±0.03 

-0.3 

1484.2960 

1486.0920 

0.81±0.02 

0.81 ±0.02 

-0.0 

0.82±0.02 

0.85±0.02 

1.2 

1501.3752 

1502.8315 

0.29±0.03 

0.25±0.03 

-1.1 

0.29±0.02 

0.26±0.02 

-0.8 

1503.6245 

1504.5574 

0.22±0.02 

0.23±0.02 

0.4 

0.24±0.01 

0.20±0.02 

-1.7 

1508.1309 

1511.7717 

0.37±0.04 

0.39±0.05 

0.3 

0.40±0.03 

0.38±0.04 

-0.3 

1513.8568 

1515.2911 

0.26±0.03 

0.20±0.03 

-1.5 

0.23±0.02 

0.23±0.02 

0.2 

1515.5824 

1517.8801 

1.47±0.02 

1.44±0.03 

-1.0 

1.46±0.02 

1.49±0.02 

1.1 

1517.8801 

1519.3525 

0.99±0.02 

0.99±0.02 

-0.2 

1.00±0.01 

1.01±0.02 

0.3 

1522.9287 

1524.4659 

0.42±0.03 

0.44±0.03 

0.4 

0.42±0.02 

0.38±0.03 

-1.2 

1525.8251 

1527.4271 

0.64±0.03 

0.62±0.03 

-0.6 

0.62±0.02 

0.63±0.03 

0.4 

1530.4724 

1531.2815 

0.21±0.02 

0.15±0.03 

-1.8 

0.18±0.02 

0.18±0.02 

-0.3 

1532.8252 

1534.3689 

0.36±0.03 

0.33±0.03 

-0.7 

0.36±0.02 

0.38±0.03 

0.4 

1536.3785 

1538.2524 

0.70±0.03 

0.71 ±0.03 

0.3 

0.73±0.02 

0.75±0.03 

0.5 

1547.2039 

1548.5242 

0.38±0.03 

0.43±0.03 

1.1 

0.38±0.02 

0.40±0.03 

0.4 

1548.8058 

1551.0970 

0.80±0.04 

0.73±0.04 

-1.1 

0.74±0.03 

0.69±0.04 

-1.1 

1551.1261 

1553.3009 

0.78±0.04 

0.72±0.05 

-1.0 

0.74±0.03 

0.72±0.04 

O 

1 

1561.0355 

1564.2556 

1.21±0.05 

1.23±0.06 

0.3 

1.18±0.04 

1.43±0.04 

4.3*’ 

1564.8542 

1565.9384 

0.31±0.03 

0.28±0.04 

-0.5 

0.29±0.02 

0.30±0.03 

0.3 

1580.8251 

1584.6600 

1.48±0.06 

1.27 ±0.09 

-1.9 

1.39±0.05 

1.69±0.05 

4.0^= 

1633.3333 

1635.2345 

0.65±0.04 

0.66±0.05 

0.1 

0.69±0.03 

0.63±0.04 

-1.1 

1635.5581 

1639.5630 

1.18±0.07 

1.07 ±0.08 

-1.1 

1.14±0.05 

I .05±0.06 

-1.1 

1644.7046 

1647.4554 

0.33±0.06 

0.39±0.07 

0.6 

0.37±0.05 

0.30±0.06 

-1.0 

1652.1115 

1655.0240 

0.50±0.06 

0.62±0.07 

1.3 

0.55±0.05 

0.56±0.06 

0.1 

1669.7450 

1672.0913 

0.90±0.05 

0.83±0.05 

-1.0 

0.85±0.04 

0.83±0.05 

-0.3 

1693.5436 

1694.5145 

0.23±0.04 

0.22±0.04 

-0.1 

0.24±0.03 

0.20±0.03 

-0.8 

1710.5782 

1713.4303 

1.03±0.07 

1.07 ±0.07 

0.5 

1.05±0.05 

1.14±0.06 

1.2 

1718.1025 

1720.0240 

0.61±0.06 

0.56±0.07 

-0.5 

0.60±0.05 

0.61±0.06 

0.2 

1726.7717 

1728.6649 

0.69±0.06 

0.71 ±0.06 

0.3 

0.67±0.05 

0.72±0.05 

0.7 

1734.0371 

1736.1002 

0.90±0.06 

0.98±0.06 

0.9 

0.92±0.05 

1.00±0.05 

1.2 

1739.0291 

1743.6892 

2.41±0.13 

2.56±0.09 

0.9 

2.47±0.07 

2.41±0.08 

-0.6 

1750.1213 

1751.5032 

0.65±0.05 

0.70±0.05 

0.7 

0.69±0.04 

0.62±0.05 

-1.1 


Notes-List of 39 strong absorption lines that are detected in Epoch-1, Epoch-2, and Epoch-3 data with EW/ctew > 5. Columns 1 and 2 are the lower and 
upper bounds of wavelengths, respectively, within which EWs are measured. Columns 3 and 4 are the observed EWs in Epoch-1 and Epoch-2 data, 
respectively. Columns 6 and 7 are the observed EWs in Epoch-1' and Epoch-3 data, respectively. The Icr errors in EW measurements (i.e. ctew) 
calculated after incorporating the continuum placement uncertainties. Columns 5 and 8 are the Sij values as defined in Eq. 1. “N IV A765, ^Ne VIII A770, 
and ^Ne VIII A780 absorption lines from the mini-BAL outflow studied here. 





